Determination of activation energy of ion-implanted deuterium release from W-Y2O3
Wang Xue-Feng1, 2, Wu Ji-Liang2, Li Qiang2, Yang Rui-Zhu2, Wang Zhan-Lei2, Chen Chang-An2, Feng Chun-Rong1, 2, Rao Yong-Chu2, Chen Xiao-Hong1, †, Ye Xiao-Qiu2, ‡
School of Sciences and Research Center for Advanced Computation, Xihua University, Chengdu 610039, China
Science and Technology on Surface Physics and Chemistry Laboratory, Mianyang 621907, China

 

† Corresponding author. E-mail: shengxiaohongb@163.com xiaoqiugood@sina.com

Project supported by the National Magnetic Confinement Fusion Energy Research Project, Ministry of Science and Technology of China (Grant No. 2015GB109002), the Innovation Fund of Postgraduate, Xihua University, China (Grant No. ycjj2018017), and the National Natural Science Foundation of China (Grant No. 21401173).

Abstract

The retention and release of deuterium in W–2%Y2O3 composite materials and commercially pure tungsten after they have been implanted by deuterium plasma (flux ∼ 3.71 × 1021 D/m2⋅s, energy ∼ 25 eV, and fluence up to 1.3 × 1026 D/m2) are studied. The results show that the total amount of deuterium released from W–2%Y2O3 is 5.23 × 1020 D/m2(2.5 K/min), about 2.5 times higher than that from the pure tungsten. Thermal desorption spectra (TDS) at different heating rates (2.5 K/min–20 K/min) reveal that both W and W–2%Y2O3 have two main deuterium trapped sites. For the low temperature trap, the deuterium desorption activation energy is 0.85 eV (grain boundary) in W, while for high temperature trap, the desorption activation energy is 1.57 eV (vacancy) in W and 1.73 eV (vacancy) in W–2%Y2O3.

1. Introduction

Tungsten has been considered as one of the candidate PFMs for ITER because of its favorable physical and chemical properties, such as high thermal conductivity, high melting point, high-temperature mechanics performance, low physical and chemical sputtering yields, and no chemical reaction with hydrogen.[17] However, the brittleness of W material leads PFM to likely undergo crack failure when the transient event occurs.[8]

W–Y2O3 alloys can effectively impede the recovery and recrystallization at high temperature, and have low ductile–brittle transition temperature (DBTT) compared with pure tungsten.[12,13] Tan et al.[14] prepared and studied the mechanical properties and microstructure changes of W–Y2O3 alloy under helium irradiation, the results showed that W–Y2O3 have the trend of lattice distortion, polycrystal and phase transformation. Yao et al.[15] studied the surface damage of W–Y2O3 irradiated by helium. The results showed that the grain orientation of W–Y2O3 irradiated by helium ion changes obviously and helium bubbles gather near the surface phase interface. The surface damage and hydrogen isotope retention of W–Y2O3 exposed to deuterium plasma have been investigated, the experimental results of Tan[16] showed that the retention of deuterium is related to the texture and laminated microstructure of the irradiated surface of the sample. The results of Zhao[17] showed that the blisters and deuterium retention of W–Y2O3 are strongly dependent on the injection temperature. However, the desorption activation energy of deuterium in W–Y2O3 is rarely reported.

In this work, based on the research of Yao et al., the deuterium desorption activation energy in W–2%Y2O3 composite material and commercial pure tungsten after being exposed to deuterium plasma are experimentally determined by thermal desorption spectroscopy. The differences in deuterium retention and release between W–2%Y2O3 and W are discussed.

2. Experimental details
2.1. Sample preparation

W–2%vol Y2O3 was supplied by Hefei University of Technology, which was prepared by using the wet chemistry combined with plasma sintering,[18,19] and the pure-W was obtained from ATTL www.tlwm.cn) with a purity of 99.95 wt%. All samples were hot rolled with the rolling ratio of 50%, and then were fabricated into square pieces with 10 mm×10 mm×1 mm from the same manufacturing batch. The surfaces were mechanically polished to a mirror-like finish then ultrasonically cleaned in acetone and ethanol. Finally, all samples were annealed at 1273 K under vacuum better than 5 × 10−5 Pa for 1 h to relieve residual stresses generated in the grinding and polishing process.

2.2. Methods

Deuterium ion irradiation was conducted by using comprehensive ECR plasma for tritium (CEPT) generating device in the Science and Technology on Surface Physics and Chemistry Laboratory.

The ion implantation direction is perpendicular to the rolling direction (RD)–transverse direction (TD) surface of the material (as shown in Fig. 1).

Fig. 1. (a) RD, TD, and ND represent rolling direction, transverse direction, and normal direction, respectively, and (b) schematic diagram of plasma implantation direction.

The detailed conditions of the deuterium ion irradiation are listed in Table 1. The pressure of deuterium (purity: 99.999%) was fixed at 0.7 Pa in the irradiation chamber during its operation. The energy of the deuterium ions and the flux were obtained by the langmuir probe. During deuterium ion irradiation, the samples were 375 K according to the measurements obtained by S-type thermocouple contacting the back of the samples.

Table 1.

Detailed experiment conditions for pure-W and W–Y2O3.

.

After the samples are irradiated by deuterium ions, the deuterium retention and desorption activation energy in the sample were obtained by using thermal desorption spectroscopy device[20] in the Science and Technology on Surface Physics and Chemistry Laboratory. The pipe which the ions passed through was baked for 72 h to remove residual gas in the material prior to the TDS experiment. The samples were heated by linearly increasing the temperature up to 1273 K at different heating rates (Table 2) in the quartz tube under vacuum better than 10−6 Pa. In addition, in order to quantitatively calculate the deuterium retention, D2 (m/z = 4) and HD (m/z = 3) signal were tracked by quadrupole mass spectrometer (QMS), and the D2/H2 signal was calibrated by using two standard leaks after experiment. The HD signal was the average of the sum of D2 and H2 signals. The x-ray diffraction (XRD) and scanning electron macroscopy (SEM) were also used to characterize the non-irradiated and irradiated samples.

Table 2.

Heating rate for pure-W and W–Y2O3 in TDS experiment.

.
3. Results and discussion
3.1. Sample characterization

Figures 2(a) and 2(b) show the orientation image maps of the W and W–Y2O3, characterized by the electron backscattered diffraction (EBSD). It can be seen from Fig. 2(c) that red color and blue color are dominant in Fig. 2(a), corresponding to (100) and (111) crystal plane, respectively; dominant in Fig. 2(b) are red color, blue color, and green color, corresponding to (100), (111), and (110) crystal planes. Figures 2(a) and 2(b) show that there are (100) and (111) textures on the RD–D surface of W and (100), (111), and (110) textures on the RD–TD surface of W–Y2O3.

Fig. 2. EBSD images of (a) W and (b) W–Y2O3 before being irradiated by deuterium ions, (c) legend for panels (a) and (b), SEM images of (d) W and (e) W–Y2O3 after being irradiated by deuterium ions.

The surface morphology of W and W–Y2O3 irradiated by deuterium ions confirm this point (Figs. 2(d)2(e))). After irradiation, the surface of W is covered with a large number of small blisters (∼ 100 nm in size, see Fig. 5(d)), The blisters on the surface of W–Y2O3 are much larger than those on the surface of W (∼ 2 μm in size), accompanied by the bulges of the substrate (Fig. 2(e)).

Fig. 3. XRD patterns of pure-W and W–Y2O3 before and after being irradiated, for (a) non-irradiated pure-W and W–Y2O3 samples; (b) pure-W before and after being irradiated; (c) W–Y2O3 before and after being irradiated; (d) irradiated pure-W and W–Y2O3 samples.
Fig. 4. TDS spectra of (a) W and (b) W–Y2O3 at different heating rates.
Fig. 5. Plots of desorption activation energies of LTT and HTT of W and W–Y2O3 by least squares fitting, where line is least square fitting.

Pan et al.[21] showed that the potential barrier for H ions to diffuse into the surface is larger than that into the bulk phase, and H ions preferentially accumulate on the W (111) surface compared with on the bulk phase, and W (110) is the resistance surface for the formation of H blisters. Figures 2(a) and 2(b) show that there are a large number of other oriented crystal planes in the red region of the RD–TD plane of W, which leads the red region to fragment, while the color boundary of the RD–TD plane of W–Y2O3 is obvious, the other orientations of the RD–TD plane are few, and the average grain boundary size is much larger than that of W. Figures 2(d) and 2(e) show that the deuterium blisters with small size (∼ 100 nm) are produced on the surface of W during irradiation, and deuterium blisters with larger size (∼ 2 μm) appear in some areas of W–Y2O3. This is consistent with Pan’s conclusion.

Figure 3 shows the XRD spectra of the main characteristic peaks of W and W–Y2O3 before and after being irradiated by the deuterium ions. In Fig. 3(a) the normal XRD peak of W in the original state, as a reference standard XRD spectrum, reflects the isotropic characteristics. As shown in Fig. 3(b), the peak position (2θ) of the irradiated W has a slight leftward shift of 0.08° compared with the peak of the non-irradiated sample, while the peak position of W–Y2O3 shifted 0.2° to the left after being irradiated by the deuterium ions. Table 3 shows the the lattice constant of W is increased by 0.009 Å while that of W–Y2O3 by 0.014 Å.

Table 3.

Lattice constants of pure-W and W–Y2O3 before and after being irradiated.

.
3.2. Deuterium release and retention

Figures 4(a) and 4(b) show the D2 TDS spectra of W and W–Y2O3. As expected, when the heating rate increases, the main desorption peaks of deuterium in W and W–Y2O3 extend to the high temperature section and the peak shape intensity is enhanced due to non-diffusion deuterium.[20]

The TDS spectra show the macroscopic characteristics of deuterium released from the whole sample, which is formed by coupling the different deuterium trap release peaks in the sample. The TDS spectra of W–Y2O3 and W have obvious shoulder peaks near 600 K, indicating that there is a corresponding low temperature deuterium trap in the sample. The whole release peak can be divided into the coupling of the release peaks of the low temperature traps (LTTs) and the high temperature traps (HTTs) (Figs. 4(a) and 4(b)).

The TDS peak of W and W–Y2O3 are in accordance with Gauss distribution. According to the results of the reaction kinetics in differential thermal analysis of Kissinger,[22] the desorption activation energy can be calculated from the TDS spectra with different heating rates by the following formula:

where β is the heating rate, Tp is the temperature of the peak of the TDS spectrum, Ea is the desorption activation energy of sample, and R is the universal gas constant.

Since β and Tp are known, the evoluation of desorption activation energy Ea of deuterium ion from a trapping site can be calculated from the slope of the versus (1/Tp) plot. The results are shown in Fig. 5. The fitting results have good linear relationship and satisfy the Arrehenius relationship. The results show that the Ea value of W–Y2O3 is 0.39 eV for LTT (∼ 600 K) and 1.73 eV for HTT (∼ 700 K), while the Ea value of pure W is 0.86 eV for LTT (∼ 600 K) and 1.57 eV for HTT (∼ 700 K) as shown in Fig. 5.

According to the experimental and modeling results (see Table 4), the values Ea of 0.39 eV and 0.86 eV are close to those explained as the migration of deuterium atoms between LTTs such as interstitials, dislocations, and grain/phase boundaries. The smaller value of Ea for LTT in W–Y2O3 is in acoordance with the fact that there are more phase boundaries in W–Y2O3 due to the doping effects. Obviously, the values of Ea of 1.57 eV and 1.73 eV for HTT can be attributed to deuterium atoms in vacancies. Interestingly, the value of Ea for HTT in W–Y2O3 is higher than that in pure W. Considering that the blisters on the surface of W–Y2O3 are larger than those of pure W, we think that there may be vacancy clusters in W–Y2O3.

Table 4.

Desorption activation energy of deuterium in W.

.

It is worth mentioning that some deuterium atoms even are desorbed from pure W at very high temperature (1200 K) (Fig. 4), while W–Y2O3 has no such behavior. The origin of such differences is not fully clear. It may be due to the fact that deuterium atoms gather on the subsurface of W–Y2O3, giving rise to blisters and leading to more defects. The capture of deuterium atoms by these defects prevents the deuterium atoms from diffuing to deeper depths. Qiu et al.[30] also found that the trapping of positions may prevent hydrogen isotope from diffusing and also from accumulating their concentration on the shallow substrate. However, for pure W, part of the deuterium atoms are captured in the sublayer; more deuterium atoms will diffuse into the bulk. Figure 6 shows that the deuterium content of W–Y2O3 is much higher than that of W in a depth range of 0 nm–2000 nm, which is in good agreement with Qiu’s conclusion.

Fig. 6. NRA determined deuterium depth profile (topmost 2000 nm) for W and W–Y2O3 exposed to deuterium plasma.

The results of Zhou et al.[31] showed that without considering the effect of temperature, most of the H atoms occupy the interstitial positions in the bulk phase, while only some H atoms occupy the defects related to the vacancies. Therefore, most of the deuterium atoms in W–Y2O3 are mainly concentrated at the interstitial positions near the surface layer, while most of the deuterium atoms in W are distributed at the interstitial positions in the whole bulk phase, which may be the reason why the deuterium release peak at the interstitial position in W–Y2O3 is observed in the TDS spectrum, but the deuterium release peak at the interstitial position in W cannot be observed obviously.

Figure 7 shows the deuterium retention of W and W–Y2O3, when considering the contribution of LTT and HTT respectively. In the W–Y2O3, the contribution of deuterium from LTT to the total amount zigzags upward with the increase of the heating rate, and finally exceeds the contribution of HTT to the total amount (Fig. 7(b)). However, for W (Fig. 7(c)), with the increase of heating rate, the deuterium release in HTT is always higher than that in LTT.

Fig. 7. (a) Deuterium retentions of W and W–Y2O3 with different heating rates; deuterium releases of LTT and HTT of (b) W–Y2O3; (c) W.
4. Conclusions

The deuterium retention and release in Y2O3 doped by tungsten, exposed to high-flux low-energy plasma, is studied. The total number of deuterium atoms released from W–Y2O3 is 5.23 × 1020 D/m2 (2.5 K/min), about 2.5 times higher than that from pure tungsten. The W–Y2O3 has two main deuterium trapped sites. For low temperature traps, the desorption activation energy is 0.39 eV while for high temperature traps (vacancy), the desorption activation energy is 1.73 eV. Both low temperature traps and high temperature traps make important contributions to trapping D.

Reference
[1] Tokunaga K Baldwin M J Doerner R P Noda N Kubota Y Yoshida N Sogabe T Kato T Schedler B 2005 J. Nucl. Mater. 337--339 887
[2] Davis J W Barabash V Makhankov A Plochl L Slattery K T 1998 J. Nucl. Mater. 256--263 308
[3] Alimov V K Hatano Y Tyburskapuschel B Sugiyama K Takagi I Furuta Y Dorner J Fuseder M Isobe K Yamanishi T Matsuyama M 2013 J. Nucl. Mater. 441 280
[4] Ou J Xiang N Men Z Z Zhang L Xu J C Gao W 2019 Chin. Phys. 28 125201
[5] Zhang H Wen S L Pan M Huang Z Zhao Y Liu Y Chen J M 2016 Chin. Phys. 25 056102
[6] Wang J Zhang H Cheng X L 2013 Chin. Phys. 22 085201
[7] Sang C F Dai S Y Sun J Z Bonnin X Xu Q Ding F Wang D Z 2014 Chin. Phys. 23 115201
[8] Wirtz M Linke J Pintsuk G Singheiser L Zlobinski M 2013 J. Nucl. Mater. 438 S833
[9] Yar M Wahlberg S Bergqvist H H Salem G Johnsson M Muhammed M 2011 J. Nucl. Mater. 412 227
[10] Veleva L Schaeublin R Battabyal M Plociski T Baluc N 2015 Int. J. Refract. Met. Hard Mater. 50 210
[11] Kim Y Lee K H Kim E Cheong D Hong S H 2009 Int. J. Refract. Met. Hard Mater. 27 842
[12] Uytdenhouwen I Decréton M Hirai T Linke J Pintsuk G Oost G V 2007 J. Nucl. Mater. 363--365 1099
[13] Smid I Akiba M Vieider G Plochl L 1998 J. Nucl. Mater. 258--263 160
[14] Tan X Luo L Chen H Zhu X Zan X Luo G Chen J Li P Cheng J Liu D Wu Y 2015 Sci. Rep. 5 12755
[15] Yao G Luo L Tan X Y Zan X Xu Q Zhu X Y Lu E Cao X Z Cheng J G Wu Y C 2019 Materialia 6 100268
[16] Tan Y Lian Y Feng F Chen Z Wang J Liu X Guo W Cheng L Lu G 2018 J. Nucl. Mater. 509 145
[17] Zhao M Jacob W Manhard A Gao L Balden M Von Toussaint U Zhou Z 2017 J. Nucl. Mater. 487 75
[18] Xu M Y Luo L M Lin J S Xu Y Zan X Xu Q Tokunaga K Zhu X Y Wu Y C 2018 J. Alloys Compd. 766 784
[19] Wang Y C Yao Z T Cheng X W Wu F Y Wang F C 2007 J. B. Inst. Thechnol. 27 824
[20] Wang Z Chen C A Song Y Xiang X Wang W Liu L Wang B 2018 Fusion Eng. 126 139
[21] Pan G Zhang Y Li Y Zhang C H Zhao Z Zeng Z 2017 Int. J. Mod. Phys. 28 1750090
[22] Kissinger H E 1957 Anal. Chem. 29 1702
[23] Degtyarenko N N Pisarev A A 2015 Phys. Procedia 71 30
[24] Eleveld H 1996 Hydrogen and helium in selected fusion reactor materials Ph.D. Dissertation Netherlands Delft University of Technology in Netherlands https://resolver.tudelft.nl/uuid:08244a63-c54b-4f91-809d-66a70bed6a7d
[25] Ogorodnikova O V Roth J Mayer M 2003 J. Nucl. Mater. 313--316 469
[26] Poon M Haasz A A Davis J W 2008 J. Nucl. Mater. 374 390
[27] Heinola K Ahlgren T 2010 Phys. Rev. 81 073409
[28] Johnson D F Carter E A 2010 J. Mater. Res. 25 315
[29] Ogorodnikova O V Roth J Mayer M 2008 J. Appl. Phys. 103 034902
[30] Qiu M Zhai L Cui J Fu B Li M Hou Q 2018 Chin. Phys. 27 073103
[31] Zhou H Liu Y Jin S Zhang Y Luo G N Lu G H 2010 Nucl. Fusion 50 115010